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We review an approach, developed over the past few years, to describe hadronic 
\ matter at finite density and temperature, whose underlying theoretical framework 

is the Skyrme model, an effective low energy theory rooted in large QCD. In 
this approach matter is described by various crystal structures of skyrmions, 
classical topological solitons carrying baryon number, from which conventional 
baryons appear by quantization. Chiral and scale symmetries play a crucial role 
in the dynamics as described by pion, dilaton and vector meson degrees of free- 
, dom. When compressed or heated skyrmion matter describes a rich phase dia- 

gram which has strong connections with the confinement /deconfinement phase 
transition. 
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1.1. Introduction 

An important issue at present is to understand the properties of hadronic matter 
under extreme conditions, e.g., at high temperature as in relativistic heavy-ion 
physics and/or at high density as in compact stars. The phase diagram of hadronic 
matter turns out richer than what has been predicted by perturbative Quantum 
Chromodynamics (QCD).^ Two approaches have been developed thus far to discuss 
this issue : on the one hand. Lattice QCD which deals directly with quark and gluon 
degrees of freedom, and on the other, effective field theories which are described in 
terms of hadronic fields. We shall describe in here a formalism for the second 
approach based on the topological soliton description of hadronic matter firstly 
introduced by Skyrme. ^'^ 

Lattice QCD, the main computational tool accessible to highly nonpcrturba- 
tive QCD, has provided much information on the the finite temperature transition, 
such as the value of the critical temperature, the type of equation of state, etc.^ 
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However, due to a notorious 'sign problem', lattice QCD could not be applied to 
study dense matter. Only in the last few years, it has become possible to simu- 
late QCD with small baryon density.^ Chiral symmetry is a flavor symmetry of 
QCD which plays an essential role in hadronic physics. At low temperatures and 
densities it is spontaneously broken leading to the existence of the pion. Lattice 
studies seem to imply that chiral symmetry is restored in the high temperature 
and/or high baryon density phases and that it may go hand- in- hand with the con- 
flnement/deconfincmcnt transition. The quark condensate (qq) of QCD is an order 
parameter of this symmetry and decreases to zero when the symmetry is restored. 

The Skyrme model, is an effective low energy theory rooted in large Nc QCD,^'^^ 
which wc have applied to dense and hot matter studies. ''^-'^^ The model does not have 
explicit quark and gluon degrees of freedom, and therefore one can not investigate 
the confinement /deconfinement transition directly, but we may study the chiral 
symmetry restoration transition which occurs close by. The schemes which aim 
at approaching the phase transition from the hadronic side are labelled 'bottom 
up' schemes. The main ingredient associated with chiral symmetry is the pion, 
the Goldstone boson associated with the spontaneously broken phase. The various 
patterns in which the symmetry is realized in QCD will be directly reflected in 
the in-medium properties of the pion and consequently in the properties of the 
skyrmions made of it. 

The most essential ingredients of the Skyrme model are the pions, Goldstone 
bosons associated with the spontaneous breakdown of chiral symmetry. Baryons 
arise as topological solitons of the meson Lagrangian. The pion Lagrangian can 
be realized non-linear ly as U = exp(iT ■ Tr/f^), which transforms as U — > QlUq^ji 
under the global chiral transformations SUL{Nf) x SUL{Nf); gL & SUL{Nf) and 
Qr € SUR{Nf). Hereafter, we will restrict our consideration to Nf = 2. In the case 
oi Nf = 2, the meson field tt represents three pions as 



The Lagrangian for their dynamics can be expanded in powers of the right and 
left invariant currents = Ud^W and = U^d^U, which transforms as Rfj, — > 
9LRng\ ^"^^ ~^ 9r^h9^r- The lowest order term is 



Here, f-rc — 93 MeV is the pion decay constant. 

Throughout this paper, wc take the following convention for the indices: (i) 
a, 6, ••• = 1,2,3 (Euchdean metric) for the isovector fields; (ii) ••• = 1,2,3 
(Euclidean metric) for the spatial components of normal vectors; (iii) ii,v, - ■ ■ = 
0, 1, 2, 3 (Minkowskian metric) for the space-time 4-vectors; (iv) a, /3, • • • = 0, 1, 2, 3 
(Euclidean metric) for isoscalar(0)-|- isovectors(l,2,3). 




(1.1) 



(1.2) 
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In the next order, one may find three independent terms consistent with Lorentz 
invariance, parity and G-parity as 

A = atr[L^, ^ /3tr{L^, L„}+ + -ftv{df,L^ f. (1.3) 

In his original work,^'^ Skyrme introduced only the first term to be denoted as 

1 



^s^. = ^tr[L^,L.V, (1-4) 

which it is still second order in the time derivatives. The value of the "Skyrme 
parameter" may be evaluated by using tttt data. In the Skyrme model, it is also 
determined, for example, as e = 5.45 to fit the nucleon-Delta masses, or as 
e = 4.75 to fit the axial coupling constant of nucleon. 

One may build up higher order terms with more and more phenomenological 
parameters. However, this naive derivative expansion leads to a Lagrangian which 
has an excessive symmetry; that is, it is invariant under U U\ which is not a 
genuine symmetry of QCD. To break it, we need the Wess-Zumino-Witten term.-^* 
The corresponding action can be written locally as 

Swzw = J d^xe^'-^f'^tviL^L, ■ ■ ■ L„), (1.5) 

i.e. in a five dimensional space whose boundary is the ordinary space and time. For 

Nf = 2 this action vanishes trivially, but for Nf = 3 it provides a hypothesized 
process KK 7r+7r°7r~. When the action is U{1) gauged for the pions to interact 
with photons, this term plays a nontrivial role even with two flavors. 

Chiral symmetry is explicitly broken by the quark masses, which provides the 
masses to the Goldstonc bosons. The mass term can be incorporated in the same 
way as chiral symmetry is broken in QCD; that is, 

jCm = ^^tr((C/ + f/t - 2)) - -Mtr(Al(i7 +U^- 2)), (1.6) 

where 

M=h-'). (1.7) 

We neglect the u- and d-quark mass difference. 

The approach has been generalized to more sophisticated meson Lagrangians 
which are constructed by implementing the symmetries of QCD.^^ The scale dilaton 
has been incorporated into the effective scheme to describe in hadronic language 
the scale anomaly. The vector mesons p and u with masses m^^^ ~ 780 MeV 
can be incorporated into the Lagrangian by using the hidden local symmetry (HLS) 

and guided by the matching of this framework to QCD in what is called 'vector 
manifestation' (VM).^"^ We shall discuss these generalizations, when required in the 
discussion of skyrmion matter, later on. 

The classical nature of skyrmions enables us to construct a dense system quite 
conveniently by putting more and more skyrmions into a given volume. Then, 
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skyrmions shape and arrange themselves to minimize the energy of the system. 
The ground state configuration of skyrmion matter are crystals. At low density 
it is made of well-localized single skyrmions.'^"' At a critical density, the system 
undergoes a structural phase transition to a new kind of crystal. It is made of 
'half-skyrmions' which are still well-localized but carry only half winding number. 
In the half-skyrmion phase, the system develops an additional symmetry which 
leads to a vanishing average value of ct = ^Tr{U), the normalized trace of the U 
field. In the studies of the late 80's,^^ the vanishing of this average value (cr) was 
interpreted as chiral symmetry restoration by assuming that (tr) is related to the 
QCD order parameter (qq). However, in Ref. 7, it was shown that the vanishing 
of (ct) cannot be an indication of a genuine chiral symmetry restoration, because 
the decay constant of the pion fluctuating in such a half-skyrmion matter does not 
vanish. This was interpreted as a signal of the appearance of a pseudogap phase 
similar to what happens in high Tc superconductors.^^ 

The puzzle was solved in Ref. 10 by incorporating a suitable degree of freedom, 
the dilaton field x, associated to the scale anomaly of QCD. The dilaton field takes 
over the role of the order parameter for chiral symmetry restoration. As the density 
of skyrmion matter increases, both (ct) and (x) vanish (not necessarily at the same 
critical density). The effective decay constant of the pion fluctuation vanishes only 
when (x) becomes zero. It is thus the dilaton field which provides the mechanism 
for chiral symmetry restoration. 

Contrary to lattice QCD, there are few studies on the temperature dependence 
of skyrmion matter. Skyrmion matter has been heated up to melt the crystal into a 
liquid to investigate the crystal-liquid phase transition ^^'^^ a phenomenon which is 
irrelevant to the restoration of chiral symmetry. We have studied skyrmion matter 
at finite density and temperature and have obtained the phase diagram describing 
the realization of the chiral symmetry. 

The contents of this review are as follows. Section [L2l deals with the history of 
skyrmion matter and how our work follows from previous investigations. We also 
study of the pion properties inside skyrmion matter at finite density. To confront the 
results with reality, in Sec. ll.Sl we show that the scale dilaton has to be incorporated 
and we discuss how the properties of the pion change thereafter. Section 11.41 is 
devoted to the study of the temperature dependence and the description of the 
phase diagram. In Sec. ll.Sl we incorporate vector mesons to the scheme and discuss 
the problem that arises due to the coupling of the w meson and our solution to it. 
Finally the last section is devoted to a summary of our main results and to some 
conclusions we can draw from our study. 
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1.2. Matter at finite density 
1.2.1. Skyrmion matter 

The Skyrme model describes baryons, with arbitrary baryon number, as static soli- 
ton solutions of an effective Lagrangian for pions.^-^ The model has been used to 
describe not only single baryon properties/^''^° but also has served to derive the 
nucleon-nucleon interaction, ^'^^ the pion-nucleon interaction,^^ properties of light 
nuclei and of nuclear matter. In the case of nuclear matter, most of the develop- 
ments^'*'^^'^^"^^ done in late 80's involve a crystal of skyrmions. 
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Fig. 1.1. Energy per single skyrmion as a function of the size parameter L. The solid circles show 
the results for massless pions and the open circles are those for massive pions. Note the rapid 
phase transition around L ~ 3.8 for massless pions. 

The first attempt to understand the dense skyrmion matter was made by 
Kutchera et al?^ These authors proceeded by introducing a single skyrmion into 
a spherical Wigner-Seitz cell without incorporating explicit information on the in- 
teraction. The presently considered conventional approaches were developed later. 
In them one assumes that the skyrmions form a crystal with a specific symmetry 
and then performs numerical simulations using this symmetry as a constraint. The 
first guess at this symmetry was made by Klebanov.^^ He considered a system 
where the skyrmions are located in the lattice site of a cubic crystal (CC) and have 
relative orientations in such a way that the pair of nearest neighbors attract max- 
imally. Goldhaber and Manton suggested that contrary to Klebanov's findings, 
the high density phase of skyrmion matter is to be described by a body-centered 
crystal (BCC) of half skyrmions. This suggestion was confirmed by numerical cal- 
culations.^^ Kugler and Shtrikman,"^^ using a variational method, investigated the 
ground state of the skyrmion crystal including not only the single skyrmion CC and 
half-skyrmion BCC but also the single skyrmion face-centered-cubic crystal (FCC) 
and half-skyrmion CC. In their calculation a phase transition from the single FCC 
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to half-skyrmion CC takes place and the ground state is found in the half-skyrmion 
CC configuration. CastiUejo et al. obtained similar conclusions. 

In Fig. 11.11 we show the energy per baryon E/B as a function of the FCC box 
size parameter L0. Each point in the figure denotes a minimum of the energy for the 
classical field configuration associated with the Lagrangians (|1.2p . (|1.4|) and (|1.6p 
for a given value of L. The solid circles correspond to the zero pion mass calculation 
and reproduce the results of Kugler and Shtrikman.'^^ The quantities L and E/B, 
appearing in the figure, are given in units of (e/^)~^ (~ 0.45fm with = 93 MeV 
and e = 4.75) and E/B in units of (67r2/^)/e 1160 MeV), respectively. The latter 
enable us to compare the numerical results on E/B easily with its Bogolmoln'y 
bound for the skyrmion in the chiral limit, which can be expressed as = 1 in 
this convention. 

In the chiral limit, as we squeeze the system from L = 6 to around L = 3.8, one 
sees that the skyrmion system undergoes a phase transition from the FCC single 
skyrmion configuration to the CC half-skyrmion configuration. The system reaches 
a minimum energy configuration at L = Lmin ~ 2.4 with the energy per baryon 
E/B ~ 1.038. This minimum value is close to the Bogolmoln'y bound for the energy 
associated to Eas. (|1.2p and (|1.4p . 

On the other hand, the configuration found at i > i,nin with the constrained 
symmetry may not be the genuine low energy configuration of the system for that 
given L. Note that the pressure P = dE/dV is negative, which implies that the 
system in that configuration is unstable. Some of the skyrmions may condense 
to form dense lumps in the phase leaving large empty spaces forming a stable 
inhomogeneous as seen in Fig. 11.11 for L ~ imin- Only the phase to the left 
of the minimum, L < Lmin, may be referred to as ^^homogeneous" and there the 
background field is described by a crystal configuration. 
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Fig. 1.2. (cr) as a function of the size parameter L. The notation is the same as in Fig. 11.11 

''A single FCC is a cube with a sidelength 2L, so that there are 4 single skyrmions in a volume of 
S-L"^, that is, the baryon number density is related to L as = 1/2L^. 
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The open circles are the solutions found with a nonvanishing pion mass, = 
140 MeVtl Comparing to the skyrmion system for massless pions, the energy per 
baryon is somewhat higher. Furthermore, there is no first order phase transition at 
low densities. 



• • • • • 




Fig. 1.3. Local baryon number densities at L = 3.5 and L = 2.0 with massive pions. For L = 2.0 
the system is (almost) a half-skyrmion in a CC crystal configuration. 

In Fig. II. 2( we show (cr), i.e. the space average value of cr as a function of 
L. In the chiral limit, (a) rapidly drops as the system shrinks and reaches zero at 
L ^ 3.8, where the system goes to the half-skyrmion phase. This phase transition 
was interpreted as a signal for chiral symmetry restoration. However, as we sail 
see in the next section, this is not the expected transition. In the case of massive 
pions, the transition in (cr) is soft. Its value decreases monotonically and reaches 
zero asymptotically, as the density increases. Furthermore, as we can see in Fig. 
11.31 where the local baryon number density is shown, for L = 2 (left) and L = 3.5 
(right) in the z — plane, the system becomes a half-skyrmion crystal at high 
density. 

Another scheme used to study multi-skyrmion systems is the procedure based on 
the Atiyah-Manton Ansatz."^^ In this scheme, skyrmions of baryon number TV are 
obtained by calculating the holonomy of Yang-Mills instantons of charge N. This 
Ansatz has been successful in describing few-nucleon systems. '^'^"^^ This procedure 
has been also applied to nuclear matter with the instanton solution on a four torus. 
The energy per baryon was found to be [E/B)min = 1-058 at Lmin — 2.47, which 



''Incorporating the pion mass into the problem introduces a new scale in the analysis and therefore 
we are forced to give specific values to the parameters of the chiral effective Lagrangian, the pion 
decay constant and the Skyrme parameter, a feature which we have avoided in the chiral limit. 
In order to proceed, we simply take their empirical values, that is, f-w = 93 MeV and e = 4.75. 
Although the numerical results depend on these values, their qualitative behavior will not. 
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Fig. 1.4. The energy cost to shift a single skyrmion from its stable position by an amount d in 
the direction of the z-axis. 

is comparable to the variational result of Kugler and Shtrikman.'^'' 

In Ref. 43, the Atiyah-Maton Ansatz is employed to get skyrmion matter from 
the 't Hooft's multi-instanton solution, which is modified to incorporate dynamical 
variables such as the positions and relative orientations of the single skyrmions. This 
description provides information on the dynamics of a single skyrmion in skyrmion 
matter. Shown in Fig. I1.4l is the energy change of the system when a single skyrmion 
is shifted from its FCC lattice site by an amount d in the direction of the z-axis. 
Two extreme cases are shown. In the case of a dense system (Lp = 2L = 5.0), the 
energy changes abruptly. For small d, it is almost quadratic in d. It implies that the 
dense system is in a crystal phase. On the other hand, in the case of a dilute system 
{Lp = 10.0), the system energy remains almost constant up to some large d, which 
implies that the system is in a gas (or liquid) phase. If we let all the variables vary 
freely, the system will prefer to change to a disordered or inhomogeneous phase in 
which some skyrmions will form clusters, as we have discussed before. 

1.2.2. Pions in Skyrmion matter 

The Skyrme model also provides the most convenient framework to study the pion 
properties in dense matter. The basic strategy is to take the static configuration 
Uo{x) discussed in Sec. II. 2. li as the background fields and to look into the properties 
of the pion fluctuating on top of it. This is the conventional procedure used to find 
single particle excitations when one has solitons in a fleld theory.^'* 

The fluctuating time- dependent pion fields can be incorporated on top of the 
static fields through the Ansatz^^ 




(1.8) 
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Skyrmion Approach to finite density and temperature 9 

where 

f/, =exp(^if-0(a;)//,^, (1.9) 

with (j) describing the fluctuating pions. 

When C/o(fO = 1(P-B = 0)i the expansion in power of (/)'s leads us to 

m = + ^m2a(f)0>'' + • • • , (1.10) 

which is just a Lagrangian for the self-interacting pion fields without any interactions 
with baryons. Here, we have written only the kinetic and mass terms relavant for 
further discussions. With a non trivial Uo{f) describing dense skyrmion matter, the 
Lagrangian becomes, 

jr ^ iG'^''(f)a^0,a^0, + ^mla{x)rr + ■■■ , (I-H) 

with 

G''\x)^a'dab + 7raTTb. (1.12) 

The structure of our Lagrangian is similar to that of chiral perturbation theory 
Eq. ()1.13|) of Refs. 46,47. These authors start with a Lagrangian containing all the 
degrees of freedom, including nucleon fields, and free parameters. They integrate 
out the nucleons in and out of an a priori assumed Fermi sea and in the process 
they get a Lagrangian density describing the pion in the medium. Their result 
corresponds to the above Skyrme Lagrangian except that the quadratic (current 
algebra) and the mass terms pick up a density dependence of the form 

- X (9'"' + MU^df^UU^d^U) + (1 - ^^p) Tr(C/ + C/t - 2), 

(1.13) 

where p is the density of the nuclear matter and 13^'^ and a are physical quantities 
obtained from the pion-nucleon interactions. Note that in this scheme, nuclear 
matter is assumed ab initio to be a Fermi sea devoid of the intrinsic dependence 
mentioned above. 

We proceed via a mean field approximation consisting in averaging the back- 
ground modifications G°'^{x) and (j{x) appearing in the Lagrangian which are re- 
duced to constants, (G"'') — GSab and (cr). Then, the Lagrangian can be rewritten 
as 

£(0*) = \d,^ld^ct>l + . • . , (1.14) 



where we have carried out a wavefunction renormalization, 0* = v G0a , which leads 
to a medium modified pion decay constant and mass as 



, /G, (1.15) 



ml _ {cr) 
VG 



(1.16) 
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Fig. 1.5. Estimates of /* //,r and m'^/mTT as functions of the baryon number density of skyrmion 
matter. 



In Fig. 11.51 we show the estimates of f^/fn and m*/m7r as a function of the 
density. As the density increases, /* decreases only to ~ 0.65/^ and then it remains 
constant at that value. Our result is different from what was the general believe:'^* 
the vanishing of (a) is not an indication of chiral symmetry restoration since the 
pion decay constant does not vanish. 

Note that (cr) 2 has the same slope at low densities, which leads to m'^/niTr ~ 1 
at low densities. Since at higher densities G becomes a constant, m*/m^ decreases 
like {(t)^^"^ with a factor which is greater than 1. As the density increases, higher 
order terms in p come to play important roles and m^jniT; decreases. A more 
rigorous derivation of these quantities can be obtained using perturbation theory.^ 

The slope of (cr) at low density is approximately 1/3. If we expand (cr) about 
/9 = and compare it with Eq. (|1.13p , we obtain 



1 - 



3 Po 



1 - 



(1.17) 



which yields SttTV ~ "2l/^/(3po) ^ 42 MeV, which is comparable with the experi- 
mental value 45 MeVH. This comparison is fully justified from the point of view of 
the expansion since both approaches should produce the same result to leading 
order in this expansion. The liner term is C(l). 

The length scale is strongly dependent on our choice of the parameters f-^ and e. 
Thus one should be aware that the p scale in Fig. 11.51 could change quantitatively 
considerably if one chooses another parameter set, however the qualitative behavior 
will remain unchanged. 



'^While this value is widely quoted, there is a considerable controversy on the precise value of this 
sigma term. In fact it can even be considerably higher than this. See Ref. 48 for a more recent 
discussion. 
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Note that the density dependence of the background is taken into account to 
ah orders. No low-density approximation, whose vahdity is in doubt except at very 
low density, is ever made in the calculation. The power of our approach is that the 
dynamics of the background and its excitations can be treated in a unified way on 
the same footing with a single Lagrangian. 

1.3. Implementing scale invariance 
1.3.1. Dilaton dynamics 

The dynamics introduced in Sec. 11.11 as an effective theory for the hadronic inter- 
actions is probably incomplete. In fact, it is not clear that the intrinsic density 
dependence required by the matching to QCD is fully implemented in the model. 
One puzzling feature is that the Wigner phase represented by the half-skyrmion 
matter with (a) = supports a non- vanishing pion decay constant. This may be 
interpreted as a possible signal for a pseudogap phase. However, at some point, 
the chiral symmetry should be restored and there the pion decay constant should 
vanish. 

This difficulty can be circumvented in our framework by incorporating in the 
standard skyrmion dynamics the trace anomaly of QCD in an effective manner.^'' 
The end result is the skyrmion Lagrangian introduced by Ellis and Lanik^^ and 
employed by Brown and Rho^^ for nuclear physics which contains an additional 
scalar field, the so called scale dilaton. 

The classical QCD action of scale dimension 4 in the chiral limit is invariant 
under the scale transformation 



under which the quark field and the gluon fields transform with the scale dimension 
3/2 and 1, respectively. The quark mass term of scale dimension 3 breaks scale 
invariance. At the quantum level, scale invariance is also broken by dimensional 
transmutation even for massless quarks, as signaled by the non- vanishing of the trace 
of the energy-momentum tensor. Equivalently, this phenomenon can be formulated 
by the non- vanishing divergence of the dilatation current -D^, the so called trace 
anomaly. 



X 



A 



X = A"\t, a > 0, 



(1.18) 




9 



(1.19) 



where P{g) is the beta function of QCD. 

Broken scale invariance can be implemented into large Nc physics by modifying 
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the standard skyrmion Lagrangian, introduced in Sec. 11.11 to 




(^] Tr(t/ + C/t-2) 



1 1 r / 1 \ 1 



(1.20) 



We have denoted the non vanishing vacuum expectation value of % as /t^,, a constant 
which describes the decay of the scalar into pions. The second term of the trace 
anomaly (jl.lOp can be reproduced by the potential energy V{x), which is adjusted 
in the Lagrangian (jl.20p so that V = dV/dx — and d^V/dx^ = at x = /x-^° 

The vacuum state of the Lagrangian at zero baryon number density is defined 
hy U — I and x = fx- The fluctuations of the pion and the scalar fields about this 
vacuum, defined through 



give physical meaning to the model parameters: as the pion decay constant, to^ 
as the pion mass, as the scalar decay constant, and as the scalar mass. For 
the pions, we use their empirical values as = 93MeV and to^ = 140MeV. We fix 
the Skyrme parameter e to 4.75 from the axial- vector coupling constant gA as in Ref. 
50. However, for the scalar field x^ no experimental values for the corresponding 
parameters are available. 

In Ref. 51, the scalar field is incorporated into a relativistic hadronic model for 
nuclear matter not only to account for the anomalous scaling behavior but also to 
provide the mid-range nuclcon-nuclcon attraction. Then, the parameters and 
are adjusted so that the model fits finite nuclei. One of the parameter sets 
is = 550 MeV and = 240 MeV (Set A). On the other hand. Song et al.^"^ 
obtain the "best" values for the parameters of the effective chiral Lagrangian with 
the "soft" scalar fields so that the results are consistent with "Brown-Rho" scaling,'*^ 
explicitly, rriy, = 720 MeV and = 240 MeV (Set B). For completeness, we consider 
also a parameter set of rriy^ = 1 GeV and ~ 240 MeV (Set C) corresponding to 
a mass scale comparable to that of chiral symmetry ^ 47r/^. 

1.3.2. Dynamics of the single skyrmion 

The procedure one has to follow can be found in Ref. 10 and is similar to the one 
discussed in Sec. 11.2.11 The first step is to find the solution for the single skyrmion 
which includes the dilaton dynamics. The skyrmion with the baryon number B — 1 
can be found by generalizing the spherical hedgehog Ansatz of the original Skyrme 
model as 



U = exp(if • (?!)//^), and x = fx + X 



(1.21) 



C/o(f) = cxp{if ■ fF{r)), and xoir) = /xC(''), 



(1.22) 
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with two radial functions F{r) and C{r). Minimization of the mass equation leads 
to a coupled set of equations of motion for these functions. In order for the solution 
to carry a baryon number, Uq has the value —1 at the origin, that is, F{x = 0) = vr, 
while there is no such topological constraint for C{x — 0). All that is required is 
that it be a positive number below 1. At infinity, the fields Uo{r) and xo(^ should 
reach their vacuum values. 




0.00 I — ' — ' — ' — ' — I — ' — ' — ' — ' — I — ' — ' — ' — ' — I — ' • ^ I II — I 
0.0 0.6 1.0 1.5 2.0 2.5 



Fig. 1.6. Profile functions F{x) and C{x) as a function of x. 

Shown in Fig. 11.61 are profile functions as a function of x{— ef-^r). F{r) and 
consequently the root mean square radius of the baryon charge show little depen- 
dence on m^. On the other hand, the changes in C{r) and the skyrmion mass are 
recognizable. Inside the skyrmion, especially at the center, C(r) deviates from its 
vaccum value 1. Note that this changes in C(r) is multiplied by in the current 
algebra term of the Lagrangian. Thus, C{r) < 1 reduces the effective f-^ inside the 
single skyrmion, which implies a partial restoration of the chiral symmetry there. 
The reduction in the effective pion decay constant is reflected in the single skyrmion 
mass. 

The larger the scalar mass is, the smaller its coupling to the pionic field and 
the less its effect on the single skyrmion. In the limit of oo, the scalar field 

is completely decoupled from the pions and the model returns back to the original 
one, where C(r) = 1, M^k = 1479 MeV and (r2)V2 ^ o.43 fm. 

1.3.3. Dense skyrmion matter and chiral symmetry restoration 

The second step is to construct a crystal configuration made up of skyrmions with 
a minimal energy for a given density. 

Referring to Refs. 10,11 for the full details, we emphasize here the role the 
dilaton field in the phase transition scenario for skyrmion matter. Let the dilaton 
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field x(r) be a constant throughout the whole space as 

x/fx = X. 



(1.23) 



Then the energy per baryon number of the system for a given density can be cal- 
culated and conveniently expressed as 

E/B{X, L) = X^{E^/B) + [EJB) + X^'^E.JB) + (2^3) [X^lnX ^ \) + \) , 

(1.24) 

where i?2, and E^ are, respectively, the contributions from the current algebra 
term, the Skyrme term and the pion mass term of the Lagrangian to the energy of 
the skyrmion system, described in Sec. II. H and (2L'^) is the volume occupied by a 
single skyrmion 

The quantity E / B{X, L) can be taken as an in medium effective potential for 
X, modified by the coupling of the scalar to the background matter. Using the 
parameter values of Ref. 11 for the Skyrme model without the scalar field, the 
effective potential E/B{X) for a few values of L behaves as shown in Fig. I1.7f a). 
At low density (large L), the minimum of the effective potential is located close 
to X = 1. As the density increases, the quadratic term in the effective potential 
E I B{X) develops another minimum &t X — Q which is an unstable extremum of the 
potential V{X) in free space. At L ^ 1 fm, the newly developed minimum competes 
with the one near X ^ 1. At higher density, the minimum shifts to X — Q where 
the system stabilizes. 





(a) 



Htm) 

(b) 



Fig. 1.7. (a) Energy per single skyrmion as a function of the scalar field X for a given L. The 
results are obtained with the {E2/B), {E4/B), and (Em/B) of Ref. 11 and with the parameter 
sets B, (b) Energy per single skyrmion as a function of L. 



In Fig. I1.7r b). we plot E / B{Xmin, L) as a function of L, which is obtained 
by minimizing E/B{X,L) with respect to X for each L. The figure in the small 
box is the corresponding value of Xmin as function of L. There we see the explicit 
manifestation of a first-order phase transition. Although the present discussion is 
based on a simplified analysis, it essentially encodes the same physics as in the more 
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rigorous treatment of x given in Ref. 10 
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Fig. 1.8. Average values of cr = ^Tr(t/) and xl fx '^^ '^^'^ lowest energy crystal configuration at a 
given baryon number density. 



We show in Fig. 11.81 the average values (ct) and {x/ fx) space for the mini- 
mum energy crystal configurations obtained by the complete numerical calculation 
without any approximation for x- These data show that a 'structural' phase transi- 
tion takes place, characterized by (cr) = 0, at lower density then the genuine chiral 
phase transition which occurs when (x) = 0. The value of (cr) becomes when 
the structure of the skyrmion crystal undergoes a change from the single skyrmion 
FCC to the half-skyrmion CC. Thus, the pseudogap phase persists in an interme- 
diate density region, where the {x/ fx) does not vanish while (cr) does.^'^ A similar 
pseudogap structure has been also proposed in hot QCD.^^ 

The two step phase transition is schematically illustrated in Fig. 11.81 Let pp 
and pc be the density at which (a) and (x) vanish, repectively. 



(A) At low density (p < pp), matter slightly reduces the vacuum value of the dilaton 
field from that of the baryon free vacuum. This implies a shrinking of the radius 
of the chiral circle by the same ratio. Since the skyrmion takes all the values 
on the chiral circle, the expectation value of a is not located on the circle but 
inside the circle. Skyrmion matter at this density is in the chiral symmetry 
broken phase. 

(B) At some intermediate densities {pp < p < pc), the expectation value of cr 
vanishes while that of the dilaton field is still nonzero. The skyrmion crystal 
is in a CC configuration made of half skyrmions localized at the points where 
cr = ±1. Since the average value of the dilaton field does not vanish, the radius 
of the chiral circle is still finite. Here, {a) = does not mean that chiral 
symmetry is completely restored. We interpret this as a pseudogap phase. 

(C) At higher density {p > pc), the phase characterized by {xl fx) — ^ becomes 
energetically favorable. Then, the chiral circle, describing the fluctuating pion 
dynamics, shrinks to a point. 
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The density range for the ocurrence of a pseudogap phase strongly depends on the 
parameter choice of nix- For small nix below 700 MeV, the pseudogap has almost 
zero size. 

In the case of massive pions, the chiral circle is tilted by the explicit (mass) 
symmetry breaking term. Thus, the exact half-skyrmion CC, which requires a 
symmetric solution for points with value a = +1 and those with ct = — 1 cannot 
be constructed and consequently the phase characterized by (a) = does not exist 
for any density. Thus no pseudogap phase arises. However, (a) is always inside the 
chiral circle and its value drops much faster than that of (x/fx)- Therefore, only if 
the pion mass is small a pseudogap phase can appear in the model. 

1.3.4. Pions in a dense medium with dilaton dynamics 

Since we have achieved, via dilaton dynamics, a reasonable scenario for chiral sym- 
metry restoration, it is time to revisit the properties of pions in a dense medium. As 
was explained in Sec. I1.2.2l and in Ref. 7, we proceed to incorporate the fluctuations 
on top of the static skyrmion crystal. (We refer to Refs. 10,11 for details.) 

Using a mean field approximation we calculate the in-medium pion mass m* 
and decay constant /* obtaining. 




, (1.25) 



ml'Zl ^( ^] a[x) ml . (1.26) 



The wave function renormalization constant gives the ratio of the in-medium 
pion decay constant /* to the free one, and the above expression arises from the 
current algebra term in the Lagrangian. The explicit calculation of to* is given in 
Ref. 11. 

In Fig. 11.91 we show the (exact) ratios of the in-medium parameters relative 
to their free-space values. Only the results obtained with the parameter set B are 
shown. The parameter set A yields similar results while set C shows a two step 
structure with an intermediate pseudogap phase. Not only the average value of xo 
over the space but also Xo(r) itself vanishes at any point in space. This is the reason 
for the vanishing of to* and /*. That is, /* really vanishes when p < pc in the 
Skyrme model with dilaton dynamics. 

At low matter density, the ratio /* / /^r can be fitted to a linear function 

^ ^ 1 - 0.24(p/po) + • • • (1.27) 

At p = po, this yields f^/ f-n = 0.76, which is comparable to the other predictions. 

In Ref. 11, the in medium modification of the x decay into two pions is also 
studied using the mean field approximation. Gathering the terms with a fluctuating 
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(a) (b) 

Fig. 1.9. (a) The ratios of the in-medium parameters to the free space parameters. The graph in a 
small box shows the masses of the pion and the scalar, (b) the in-medium decay width F* (x tttt) 
as a function of p. 



scalar field and two fluctuating pion fields, we get the Lagrangian density for the 
process x ^ tttt 

^M,x^2 = ^{6ab + gab)xd,,(l)ad''(j)b, (1.28) 

where only the term from Ca is used. Averaging the space dependence of 
the background field configuration modifies the coupling constant by a factor 
((Xo//x)(l + = ((Xo//x)(l - fi'^))- Taking into account the appropriate 

wave function renormalization factors, Z^r, and the change in the scalar mass, we 
obtain the in-medium decay width as 



r*(x 



3m 



327r/. 



((Xo//x)(l- V)> 



3m 



327r/. 



*2 ■ 



(1.29) 



((Xo//x)^(l-|^^)) 

We show in Fig. I1.9l the in-medium decay width predicted with the parameter set 
B. In the region p > ppt where xo = 0, F* cannot be defined to this order. Near the 
critical point, the scalar becomes an extremely narrow-width excitation, a feature 
which has been discussed in the literature as a signal for chiral restoration. ^^'^^ 

Another interesting change in the properties of the pion in the medium is asso- 
ciated with the in medium pion dispersion relation. This relation requires, besides 
the mass, the so-called in medium pion velocity, v^^. This property allows us to 
gain more insight into the real time properties of the system under extreme con- 
ditions and enables us to analyze how the phase transition from normal matter 
to deconfined QCD takes place from the hadronic side , the so called 'bottom up' 
approach. 

At nonzero temperature and/or density, the Lorentz symmetry is broken by the 
medium. In the dispersion relation for the pion modes (in the chiral limit) 

pI^vUpI^ (1.30) 
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the velocity v^^ which is 1 in free-space must depart from 1. This may be studied 
reUably, at least at low temperatures and at low densities, via chiral perturbation 
theory.^^ The in-medium pion velocity can be expressed in terms of the time com- 
ponent of the pion decay constant, fl and the space component, /^,^^'^^ 



(0|Ai|7r''(p))i„.,„edium = */^<5«V, 



bf^\\ ^Afsxab^i (1-31) 



The conservation of the axial vector current leads to the dispersion relation (|1.30p 
with the pion velocity given by 

vl = f^/n- (1.32) 

In Ref. 60 two decay constants, ft and fs, are defined differently from those of 
Ea. (|1.3ip ). through the effective Lagrangian, 

£eff = ^TridoU^doU) - ^Tr(a,C/t9,{7) + • • • , (1.33) 

where U is an SU(2)-valued chiral field whose phase describes the in-medium pion. 
In terms of these constants, the pion velocity is defined by 

V. = fs/ft- (1.34) 

In Ref. 11, it is shown that local interactions with background skyrmion matter 
lead to a breakdown of Lorcntz symmetry in the dense medium and to an effective 
Lagrangian for pion dynamics in the form of Eq. (|1.33|) . The results are shown in 
Figs. 11.101 Both of the pion decay constants change significantly as a function of 
density and vanish - in the chiral limit - when chiral symmetry is restored. However, 
the second-order contributions to the fs and which break Lorentz symmetry, 
turn out to be rather small, and thus their ratio, the pion velocity, stays ~ 1. 
The lowest value found is ~ 0.9. Note, however, the drastic change in its behavior 
at two different densities. At the lower density, where skyrmion matter is in the 
chiral symmetry broken phase, the pion velocity decreases and has the minimum 
at p = Pp. If one worked only at low density in a perturbative scheme, one would 
conclude that the pion velocity decreases all the way to zero. However, the presence 
of the pseudogap phase transition changes this behavior. In the pseudogap phase, 
the pion velocity not only stops decreasing but starts increasing with increasing 
density. In the chiral symmetry restored phase both ft and vanish. Thus their 
ratio makes no sense. 

In Refs. 12,13, the in-medium modification of the neutral pion decay process into 
two gammas and neutrino-anti-neutrino pair are studied in the same manner. The 
tt" —^ 77 process is shown to be strongly suppressed in dense medium, while the 
process tt^ — > lyu forbidden in free space becomes possible by the Lorentz symmetry 
breaking effect of the medium. 
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Fig. 1.10. In-medium pion decay constants and their ratio, the pion velocity. 

1.4. Skyrmion matter at finite temperature 

There are many studies of lattice QCD at finite temperature. The situation is com- 
pletely different for skyrmion matter where the number of studies is limited. For 
example, skyrmion matter has been heated up to melt the crystal into a liquid to 
study the crystal- liquid phase transition. ^^--^^ However this phenomenon is irrele- 
vant to study the restoration of chiral symmetry, which interests us here for the 
reasons discussed in previos sections. 

What happens if we heat up the system? Naively, as the temperature increases, 
the kinetic energy of the skyrmions increases and the skyrmion crystal begins to 
melt. The kinetic energy associated with the translations, vibrations and rotations 
of the skyrmions is proportional to T. This mechanism leads to a solid-liquid-gas 
phase transition of the skyrmion system. However, we are interested in the chiral 
symmetry restoration transition, which is not related to the melting. Therefore, a 
new mechanism must be incorporated to describe chiral symmetry restoration. We 
show in what follows that the thermal excitation of the pions in the medium is the 
appropriate mechanism, since this phenomenon is proportional to T** and therefore 
dominates the absorption of heat. 

The pressure of non-interacting pions is given by 

P=^T\ (1.35) 

where we have taken into account the contributions from three species of pion, 
TT"*", 7r°, 7r~. This term contributes to the energy per single skyrmion volume as 
ZPV{x/ fxf' ■ The kinetic energy of the pions arises from Ca (II. 2p . and therefore 
scale symmetry implies that it should carry a factor x^. The factor 3 comes from 
the fact that our pions are massless. 
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To estimate the properties of skyrmion matter at finite temperature let us take 
X as a constant field as we did in Sec. 11.3.31 After including thermal pions, Eq. 
(|1.24p can be rewritten as 



E/Bip,T,X)= {E2/B){p) + ^T^V]X^+{Ei/B)ip)+X\\nX~l) + ^), (1.36) 



10 

where we have dropped the pion mass term. 

As in Sec. 11.3.31 chiral restoration will occur when the value of X^nin that 
minimizes E/B vanishes. By minimizing E/B with respect to X, we observe that 
the phase transition from a non-vanishing X = e^^/^ to X = 0. Thus, the nature 
of the phase transition is of the first order. 

After a straightforward calculation we obtain, 

pWi?) + ^T,^ = |9f. (1.37) 

which leads to 

Tc=[i^ [^-p'iE,/B){p^^)]] (1.38) 




For p = (zero density), our estimate for the critical temperature is 



^=-(^^1 - 205 MeV, (1.39) 



where we have used the following values for the parameters, fx = 210 MeV and 
= 720 MeV. It is remarkable that our model leads to Tc ~ 200 MeV, which is 

close to that obtained by lattice QCD ^ and in agreement with the data.^^ To us 

this is a confirmation that the mechanism chosen for the absorption of heat plays a 

fundamental role in the hadronic phase. 

The numerical results on E2/B that minimize the energy of the system for a 

given pb can be approximated by 

E./B-\l%Y'^'l'' (1.40) 
[36/^/esk, p<Po, 

where pQ = (esk/7r/3.6)^. 

Using Ea. (|1.40p for E2/B, we obtain the critical density for chiral symmetry 
restoration at zero temperature as 

P'iT - 0) = (^^Yo^ j - 0.37 fm-^ (1.41) 

Since po = 0.24 fm^'^ < p'^{T — 0) our result is consistent with the high density 
formula for E2/B used. 

The resulting critical density p''{T = 0) 0.37 fm~^ is only twice normal nuclear 
matter density and it is low with respect to the expected values. This result does 
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Fig. 1.11. The chiral phase transition. The solid line shows the exact calculation, while the gray 
lines two approximate estimates. 



not represent a problem since p'^{T — 0) scales with (/x"^x//7r) and with 
(/x'^x)^^^ and small changes in the parameters lead to larger values for the critical 
density without changing the critical temperature too much. 

For a finite density smaller than p'^iT = 0), we obtain the corresponding critical 
temperature by substituting the asymptotic formulas (|1.40p for E2/B, 



T,{p = 0) (1 - 3.09 pcf/^ for p < po, 
T,{p - 0) (1 - 1.92 p2/3)i/4 for p > 



(1.42) 



where the density is measured in fm~"^ The gray lines in Fig. 11.111 show these two 
curves. The results from the exact calculations obtaned by minimization of the 
energy (|1.36p are shown by black dots connected by black line in Fig. 11.111 The 
resulting phase diagram has the same shape but the values of the temperatures and 
densities are generally smaller than in the approximate estimates. 



1.5. Vector mesons and dense matter 



In our effort to approach the theory of the hadronic interactions and inspired by 
Weinberg's theorem we proceed to incorporate to the model the lowest-lying 
vector mesons, namely the p and the lo. In this way we also do away with the ad 
hoc Skyrme quartic term. It is known that these vector mesons play a crucial role 
in stabilizing the single nucleon system ^^'^^ as well as in the saturation of normal 
nuclear matter. 

We consider a skyrmion-type Lagrangian with vector mesons possessing hidden 
local gauge symmetry, spontaneously broken chiral symmetry and scale symme- 



try.' 



Such a theory might be considered as a better approximation to reality 
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than the extreme large approximation to QCD represented by the Skyrme model. 
Specifically, the model Lagrangian, which we investigate, is given by 

C = :!ff-^\ TT{^^U^^^'U) + ^^^(^^ Tr(C/ + C/t-2) 
4 \/x/ ^ \jx/ 

■ / \ ^ 



4 ' 



+ \goj^B^ + \d^xd'^x - [(x//x)'(ln(x//x) - i) + 3] , (1-43) 

where, U = exp(if • vr//^) = = "^^S^C- "^1^ = Pm<^ = ^^Pi. - c'^Pm + 

gp^ X p„ = d^u:, - d,u;^, and = 24^e'^''"'3Tr([/ta^[/c/ta^[/c/ta^f/). Note 

that the Skyrme quartic term is not present. The vector mesons, p and lo, are 
incorporated as dynamical gauge bosons for the local hidden gauge symmetry of 
the non-linear sigma model Lagrangian and the dilaton field x is introduced so that 
the Lagrangian has the same scaling behavior as QCD. The physical parameters 
appearing in the Lagrangian are summarized in Table 11.11 . 



Table 1.1. Parameters of the model Lagrangian 



notation 


physical meaning 


value 


u 


pion decay constant 


93 MeV 


fx 


X decay constant 


210 MeV 


9 


pTTTT coupling constant 


5.85* 


m-^ 


pion mass 


140 MeV 




X meass 


720 MeV 


mv 


vector meson masses 


770 MeVt 


a 


vector meson dominance 


2 



'v- = '"p = ™-i = o-fia with 
a = 2. cf. gpirir = 6.11 from the decay width of p — > tttt. 
t experimentally measured values are nip =768 MeV and mtj=782 MeV. 



1.5.1. Dynamics of the single skyrmion 

The spherically symmetric hedgehog Ansatz for the B = 1 soliton solution of the 
standard Skyrme model can be generalized to 

U^=^ =cxp{iT-fF{r)), (1.44) 

p;:Z^'=e^>^^r''^, p;:ur=0, (1.45) 

a;^=i = 0, c.^=oi - UWir), (1.46) 

X^=^=/xC(r). (1.47) 
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Fig. 1.12. Profile functions - F(r), G(r), W{r) and C(r). 



The boundary conditions that the profile functions satisfy at infinity are 

F{oo) = G{oo) = W(oo) = 0, C(oo) = 1, (1.48) 
and at the center (r = 0) are 

F(0) = TT, G(0) -2, VK'(O) = C"(0) ^ 0. (1.49) 

The profile functions are obtained numerically by minimizing the soliton mass 
with the boundary conditions (see Ref. 8 for the technical details). The results 
are summarized in Table 11.21 and the corresponding profile functions are given in 
Fig. 11.121 The role of the uj meson that provides a strong repulsion is prominent. 
Comparing the -np model with the -Kpuj model, the presence of the lo increases the 
mass by more than 415 MeV and the size, i.e. (r^), by more than .28 fm^. 

How does the dilaton affect this calculation? The irp model with much smaller 
skyrmion has a larger baryon density near the origin and this affects the dilaton, 
significantly changing its mean-field value from its vacuum one. The net effect of 
the dilaton mean field on the mass is a reduction of ~ 150 MeV, whereas for the 
TTpoj model it is only of 50 MeV. The details can be seen in Table [L2l The effect on 
the soliton size is, however, different: while the dilaton in the -np model produces an 
additional localization of the baryon charge and hence reduces (r^) from .21 fm^ to 
.19 fm^, in the irpuj model, on the contrary, the dilaton produces a delocalization 
and increases (r^) from .49 fm^ to .51 fm^. We will see, however, that this strong 
repulsion provided by uj causes a somewhat serious problem in the chiral restoration 
of the skyrmion matter at higher density. 
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Tabic 1.2. Single skyrmion mass and various contributions to it. 



Model 


(r2> 


£.13 = 1 




1 






eS^'' 


I?S=1 


E^=' 


7rp-model 


0.27 


1054.6 


400.2 4 


- 9.2 


110.4 


534.9 


0.0 


0.0 


0.0 


7rpx-™odel 


0.19 


906.5 


103.1 4 


- 1.4 


155.1 


504.1 


0.0 


0.0 


142.8 


7rp(i)-model 


0.49 


1469.0 


767.6 + 


39.9 


33.2 


370.7 


-257.6 


515.1 


0.0 


Trpiiix-niodel 


0.51 


1408.3 


646.0 + 


29.2 


34.9 


355.7 


-278.3 


556.7 


64.2 




Fig. 1.13. E/B, (x) and (cr) as a function of L in the models (a) without the ui and (b) with 
omega. 



1.5.2. Skyrmion Matter : an FCC skyrmion crystal 

Again, the lowest-energy configuration is obtained when one of the skyrmions is 
rotated in isospin space with respect to the other by an angle tt about an axis 
perpendicular to the line joining the two.* If we generalize this Ansatz to many- 
skyrmion matter, we obtain that the configuration at the classical level for a given 
baryon number density is an FCC crystal where the nearest neighbour skyrmions 
are arranged to have the attractive relative orientations.^ Kugler's Fourier series 
expansion method'^'' can be generalized to incorporate the vector mesons, although 
some subtelties associated with the vector fields have to be implemented. The 
details can be found in Ref. 8. 

Figs. II. 131 are the the numerical results of the energy per baryon E/B, (x) and 
(a) in various models as a function of the FCC lattice parameter L. In the npx 
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model, as the density of the system increases (L decreases), E / B changes little. It is 
close to the energy of a B = 1 skyrmion up to a density greater than po {L 1.43). 
This result is easy to interpret. As we discussed before the size of the skyrmion in 
this model is very small and therefore the skyrmions in the lattice will interact only 
at very high densities, high enough for their tails to overlap. 

In the absence of the to, the dilaton field plays a dramatic role. A skyrmion mat- 
ter undergoes an abrupt phase transition at high density at which the expectation 
value of the dilaton field vanishes (x) = 0. (In general, (x) — does not necessarily 
require (x^) — 0. However, since x ^ 0; (x) = always accompanies x = in the 
whole space.) The p meson on the other hand is basically a spectator at the classi- 
cal level, producing little change with respect to our previously studied ttx model 
except that at high densities, once the p starts to overlap, the energy of nuclear 
matter increases due to its the repulsive effect at short distances. The densities 
have to be quite high since these skyrmions are very small. Since x vanishes at the 
phase transition, we recover the standard behavior, namely, f*—0 and m* = 0. 

In the TT/awx model, the situation changes dramatically. The reason is that 
the Lu provides not only a strong repulsion among the skyrmions, but somewhat 
surprisingly, also an intermediate range attraction. Note the different mass scales 
between Figs. ll-lSf a) and ll.lSF b). In both the irpu! and the Trpux models, at high 
density, the interaction reduces E/B to 85% of the B = 1 skyrmion mass. This 
value should be compared with 94% in the np model. In the Trpx-niodel, E/B goes 
down to 74% of the B = 1 skyrmion mass, but in this case it is due to the dramatic 
behavior of the dilaton field. 

In the npujx model the role of the dilaton field is suppressed. It provides a only a 
small attraction at intermediate densities. Moreover, the phase transition towards 
its vanishing expectation value, (x) — 0, does not take place. Instead, its value 
grows at high density! 

The problem involved is associated with the Lagrangian (I1.43P which includes 
an anomalous part known as Wess-Zumino term, namely the coupling of the w to 
the Baryon current, B^. To see that this term is the one causing the problem, 
consider the energy per baryon contributed by this term.* 



where "Box" corresponds to a single FCC cell. Note that while the integral over x 
is defined in a single FCC cell, that over x' is not. Thus, unless it is screened, the 
periodic source Bq filling infinite space will produce an infinite potential w which 
leads to an infinite {E/ B)wz ■ The screening is done by the omega mass, to* . Thus 
the effective w mass cannot vanish. Our numerical results reflect this fact: at high 
density the Bq-Bq interaction becomes large compared to any other contribution. 
In order to reduce it, x has to increase, and thereby the effective screening mass 
TO* TO(^(x) becomes larger. In this way we run into a phase transition where the 
expectation value of x does not vanish and therefore /jr does not vanish but instead 




(1.50) 
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increases. 



1.5.3. A Resolution of the u) problem 

Assuming that there is nothing wrong with (|1.43p . we focus on the Wess-Zumino 
term in the Lagrangian. Our objective is to find an alternative to (jl.43p that 
leads to a behavior consistent with the expected behavior. In the absence of any 
reliable clue, we try the simplest, admittedly ad hoc, modification of the Lagrangian 
(I1.43P that allows a reasonable and appealing way-out.^ Given our ignorance as to 
how spontaneously broken scale invariance manifests in matter, we shall simply 
forego the requirement that the anomalous term be scale invariant and multiply the 
anomalous u ■ B term by {xl fx)"" n>2. We have verified that it matters little 
whether we pick n = 2 or n = 3.^ We therefore take n 3: 



L'an = %9{x/fxf^,B>^ (1-51) 
This additional factor has two virtues: 

i) It leaves meson dynamics in free space (i.e. xl fx ~ 1) unaffected, since chiral 
symmetry is realized a la sigma model as required by QCD. 

ii) It plays the role of an effective density-dependent coupling constant so that at 
high density, when scale symmetry is restored and x/ fx ~^ 0' there will be 
no coupling between the lo and the baryon density as required by hidden local 
symmetry with the vector manfiestation. 

The properties of this Lagrangian for the meson {B = 0) sector are the same as 
in our old description. The parameters of the Lagrangian are determined by meson 
physics as given in Table 11.11 

Fig. 11.141 summarizes the consequences of the modification. Depending on the 
dilaton mass, the properties of a single skyrmion show distinguished characters 
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and consequently undergoes different phase transition. A small dilaton mass, say 
m-^ 1 GeV, leads to a very small skyrmion with an rms radius about 0.1 fm. The 
light dilaton seems to react quite sensitively to the presence of the matter. One can 
see that at the center of the single skyrmion the chiral symmetry is almost restored. 
It weakens most of the repulsion from uj — B coupling, which leads us to such a small 
sized skyrmion. Since these small skyrmions are already chiral-symmetry-restored 
objects, simply filling the space with them restores the symmetry. As shown in 
the small box, chiral symmetry is restored simultaneously when (cr) vanishes. In 
case of having a large mass, the dilaton does not play any significant role in the 
structure of a single skyrmion. This scenario leads, as the density of skyrmion 
matter increases, first to a pseudogap phase transition where (cr) — and thereafter, 
at higher density, to a genuine chiral symmetry restoration phase transition where 
ix/fx) ~ ^- ^^y^^Y: whether the dilaton is light or heavy, we finally have a 
reasonable phase transition scenario that at some critical density chiral symmetry 
restoration occurs where {x/ fx) vanishes. 

Under the same mean field approximation, this skyrmion approach to the dense 
matter leads us to the scaling behaviors of the vector mesons 



With a = 1, a remarkably simple BR scaling law is obtained. These scaling laws 
imply that as the density of the matter increases the effecive quantities in medium 
scale down. 

We have shown how a slight modification of the Lagrangian resolves the lu prob- 
lem. However, in modifying the Lagrangian we have taken into account only the 
phenomenological side of the problem. Multiplying the Wess-Zumino term by the 
factor ix/fx)"' ^'^ sound theoretical support. It breaks explicitly the scale in- 
variance of the Lagrangian. Recall that the dilaton field was introduced into the 
model to respect scale symmetry. Furthermore, we don't have any special rea- 
son for choosing n = 3, except that it works well. Recently, a more fundamental 



explanation for the behavior in Eq. (|1.5ip has been found. |^ 
1.6. Conclusions 

In trying to understand what happens to hadrons under extreme conditions, it is 
necessary that the theory adopted for the description be consistent with QCD. In 
terms of effective theories this means that they should match to QCD at a scale close 

"^Private communication by M. Rho on work in progress by H.K. Lee and M. Rho. 




(1.52) 



while that of the pion decay constant is 




(1.53) 
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to the chiral scale ^ 47r/^ ~ 1 GeV. It has been shown that this matching can 
be effectuated in the framework of hidden local symmetry (HLS) and leads to what 
is called 'vector manifestation' (VM) which provides a theoretical support for a 
low-energy effective field theory for hadrons and which gives, in the chiral limit, an 
elegant and unambiguous prediction of the behavior of light-quark hadrons at high 
temperature and/or at high density. Following the indications of the HLS theory, 
we have described a Skyrmc model in which the dilaton field X) whose role in dense 
matter was first pointed out by Brown and Rho,'''^ and the vector meson fields p 
and Lo were incorporated into the Skyrme Lagrangian to construct dense skyrmion 
matter. 

We have presented an approach to hadronic physics based on Skyrme's philoso- 
phy, namely that baryons are solitons of a theory described in terms of meson fields, 
which can be justified from QCD in the large Nc expansion. We have adopted the 
basic principles of effective field theory. Given a certain energy domain we describe 
the dynamics by a Lagrangian defined in terms of the mesonic degrees of freedom 
active in that domain, we thereafter implement the symmetries of QCD and VM, 
and describe the baryonic sectors as topological winding number sectors and solve 
in these sectors the equations derived from the Lagrangian with the appropriate 
boundary conditions for the sector. In this way one can get all of Nuclear Physics 
out of a single Lagrangian. We have studied the B=l sector to obtain the proper- 
ties of the single skyrmion, the B=2 sector to understand the interaction between 
skyrmions, and our main effort has been the study skyrmion matter, as a model for 
hadronic matter, investigating its behavior at finite density and temperature and 
the description of meson properties in that dense medium. 

Skyrmc models have been proven successful in describing nuclei, the nucleon- 
nucleon interaction and pion-nucleon interactions. It turns out that Skyrme models 
also represent a nice tool for understanding low density cold hadronic matter and 
the behavior of the mesons in particular the pion inside matter. We have shown 
in here that when hadronic matter is compressed and/or heated Skyrme models 
provide useful information on the chiral phase transitions. Skyrmion matter is 
realized as a crystal and we have seen that at low densities it is an FCC crystal 
made of skyrmions. The phase transition occurs when the FCC crystal transforms 
into a half skyrmion CC one. In our study we have discovered the crucial role of the 
scale dilaton in describing the expected phase transition towards a chiral symmetry 
restored phase. We have also noticed the peculiar behavior of the iv associated to 
its direct coupling to the baryon number current and we have resolved the problem 
by naturally scaling the coupling constant using the scale dilaton. 

Another aspect of our review has been the study of the properties of elementary 
mesons in the medium, in particular those involved in the model, the pion and the 
dilaton. Moreover we have described how their properties change when we move 
from one phase to another. 

A description of the chiral restoration phase transition in the temperature- 
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density plane has been presented, whose main ingredient is that the dominant 
scenario is the absorption of heat by the fluctuating pions in the background of 
crystal skyrmion matter. This description leads to a phase transition whose dy- 
namical structure is parameter independent and whose shape resembles much the 
conventional confinement/deconfiment phase transition. We obtain, for parameter 
values close the conventional ones, the expected critical temperatures and densities. 

For clarity, the presentation has been linear, in the sense, that given the La- 
grangian we have described its phenomenology, and have made no effort to interpret 
the mechanisms involved and the results obtained from QCD. In this way we have 
taken a 'bottom up' approach: the effective theory represents confined QCD and it 
should explain the hadronic phenomenology in its domain of validity. 

The main result of our calculation is the realization that phase transition sce- 
nario is not as simple as initially thought but contains many features which make 
interesting and phenomcnologically appealing. It is now time to try to collect ideas 
based on fundamental developments and see how our effective theory and the prin- 
ciples that guide it realize these ideas. In this line of thought, it is exciting to have 
unveiled scenarios near the phase transition of unexpected interesting phenomenol- 
ogy in line with recent proposals. ^^'^^ 
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